close all
clear all
drawnow

a = linspace(0, 2^15, 128);
b = sqrt(a);

figure
plot(a, b, 'k-');

p = polyfit(a, b, 1);

hold on
plot(a, polyval(p, a), 'g--');